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Contrary to young neutron stars, young strange stars are 
not subject to the r-mode instability which slows rapidly ro- 
tating, hot neutron stars to rotation periods near 10 ms via 
gravitational wave emission. Young millisecond pulsars are 
therefore likely to be strange stars rather than neutron stars, 
or at least to contain significant quantities of quark matter in 
the interior. 
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It has recently been shown, that the emission of grav- 
itational radiation due to r-mode instabilities in hot, 
young neutron stars severely limits the rotation period 
of these stars Jj]-E|- Within the first year, where the star 
cools to 10 9 K, the star will spin down to a few per cent 
of the Kepler-limit (the ultimate rotation frequency 
where the stellar equator "falls off"), corresponding to 
a rotation period of 10-20 milliseconds. This scenario 
has several interesting ramifications First of all, 

the "slow" rotation predicted for young pulsars seems to 
agree with the few observed periods of pulsars known 
to be young, such as Crab (the rapid millisecond pul- 
sars with periods as small as 1.56 ms are expected to be 
old, cold pulsars spun-up by angular momentum accre- 
tion in binary systems, and are therefore not subject to 
the r-mode instability). Second, the spin-down involves 
emission of up to 10 52 ergs of gravitational waves, mak- 
ing the radiation potentially observable. And third, the 
instability rules out a scenario where millisecond pulsars 
form by accretion-induced collapse of a white dwarf, be- 
cause the temperature during such a process will be high 
enough for the r-mode instability to lead to spin-down. 

All rotating relativistic stars are generically unsta- 
ble against the r(otational)-mode instability [0J^]. To 
first order in the rotational frequency O of the star, 
the frequency observed by a distant inertial observer is 
LUi — u) r — mil, where the frequency in the rotating frame 
is io r — 2mfl/l(l + 1) for a mode with spherical harmonic 
indices (I, rri). Thus modes which are counterrotating in 
the comoving frame will appear corotating to a distant 
observer. Therefore, r-modes are generically unstable to 
the emission of gravitational waves; the emission removes 
positive angular momentum from a mode with nega- 
tive angular momentum in the corotating frame, thereby 
making the angular momentum and energy increasingly 
negative. In contrast to other known modes, the r-mode 
instability does not require a large rotation rate, but sets 
in even as the rotation frequency goes to zero, at least 



in a non-viscous system (for a review of instabilities in 
rotating relativistic stars the reader is referred to ||). 

Important for the r-mode instability scenario is the 
effect of internal fluid dissipation in the star that tends to 
suppress the instabilities. Lindblom, Owen, and Morsink 
||, and Andersson, Kokkotas, and Schutz Q] have shown 
that for typical neutron star equations of state, the shear 
viscosity limits the r-mode instability at temperatures 
below 10 9 K, whereas bulk viscosity is the limiting factor 
above this temperature. And they find that the rotation 
rate for standard assumptions about neutron star cooling 
is reduced to 7-8% of the Keplerian limit. 

Whereas the shear viscosity of quark matter || is 
roughly comparable to that of neutron star matter for 
the parameter range of interest here, the bulk viscosity 
of quark matter is larger by many orders of magnitude 
JlO[ | . Therefore, one should expect the r-mode instabil- 
ity to be significantly suppressed if the stars are made 
up of u-, d-, and s-quarks (so-called strange stars 
rather than ordinary neutron star matter. A similar re- 
sult would apply if quark matter is metastable rather 
than absolutely stable, so that a neutron star contains 
quark phase in its central regions. 

In the following I will show that this is indeed the case. 
The r-mode instability does not play any role in a young 
strange star. Therefore, finding a young pulsar with a ro- 
tation period below 5-10 milliseconds would be a strong 
indication of the existence of strange stars, or at least a 
significant quark content in neutron stars, and therefore 
of the stability or metastability of bulk strange quark 
matter. In addition it would re-introduce the possibility 
of forming some millisecond pulsars by accretion-induced 
collapse of white dwarf stars. Unfortunately it would also 
remove r-mode instability braking of young neutron stars 
as a strong source of gravitational radiation. 

The shear viscosity of strange quark matter due to 
quark scattering was calculated in |9). The results for 
T <C fJ,, where T is the temperature, and fi w 300MeV 
the quark chemical potential, can be written as 
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where as is the fine structure constant of the strong in- 
teractions, T g ee T/10 9 K, and p 15 = p/lO^gcm" 3 . 

The bulk viscosity of strange quark matter |l(J de- 
pends mainly on the rate of the non-leptonic weak inter- 
action [ftl 
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Because the strange quark is much more massive than 
up and down quarks, this reaction changes the concen- 
trations of down and strange quarks in response to the 
density changes involved in vibration or rotational insta- 
bilities, thereby causing dissipation. This dissipation is 
most efficient if the rate of reaction (|2|) is comparable to 
the frequency of the density change. If the weak rate is 
very small, the quark concentrations keep their original 
values in spite of a periodic density fluctuation, whereas 
a very high weak rate means that the matter immediately 
adjusts to follow the true equilibrium values reversibly. 
But in the intermediate range dissipation due to PdV- 
work is important. 

The importance of dissipation due to Eq. (||) was first 
stressed by Wang and Lu jl3) in the case of neutron stars 
with quark cores. These authors made a numerical study 
of the evolution of the vibrational energy of a neutron 
star with an O.2M quark core, governed by the energy 
dissipation due to Eq. (|J). Sawyer [[[ij expressed the 
damping in terms of the bulk viscosity, a function of tem- 
perature and oscillation frequency, which he tabulated 
for a range of densities and strange quark masses. This 
tabulation has later been used in studies of quark star 
vibration , and of the gravitational radiation reaction 
instability in rapidly rotating quark stars Jl6| . The latter 
study concluded, that the bulk viscosity was large enough 
to be important for temperatures exceeding 10 8 K, but 
that it should be a few orders of magnitude larger to 
generally dominate the stability properties. 

However, as pointed out in Ref. pp| , the bulk viscosi- 
ties in depend on the assumption, that the rate of 
Eq. (Q) can be expanded to first order in Sfi = fj, s — fid, 
where fit w 300MeV are the quark chemical poten- 
tials. This assumption is not correct at low temperatures 
(2ttT <C Sfi), where the dominating term in the rate is 
proportional to <5/x 3 . Furthermore, the rate in fbil is too 
small by an overall factor of 3, and another discrepancy 
of 2-3 orders of magnitude, perhaps due to unit conver- 
sions, appears as well. Taken together, these effects led to 
an upward correction of the bulk viscosity by several or- 
ders of magnitude, and thereby increases the importance 
for the astrophysical applications. The non-linearity of 
the rate also means, that the bulk viscosity is no longer 
independent of the amplitude of the density variations. 
The resulting bulk viscosity is (in cgs-units, with strange 

1/3 

quark mass m s , T, and [i d w 359MeVp^ in MeV, and 
the frequency of the perturbation lu in s _1 ) 
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For typical values (m s = 100 MeV, [i d = 300MeV, 
lu = lu,. = 2mfl/l(l + 1) = 2 x 10 3 s -1 ) this is 
C w 1.6 x 10 30 [93(Au/i; ) 2 + 0.29T 9 2 ] gem^s" 1 , where 
Av/vq is the perturbation amplitude. 



Kokkotas and Stergioulas || give an analytic descrip- 
tion of the r-mode instability in uniform density neutron 
stars, finding results in good agreement with more so- 
phisticated numerical calculations. The approximation 
of constant density is even better for a strange star, ex- 
cept very close to the gravitational instability limit, so 
for the estimates in the present investigation it is safe to 
use their approximate formulae for damping times etc. 
In particular, the damping time due to shear viscosity is 
given by || 

T = ? *L (4) 

sv 4tt(Z - 1)(2Z + 3) r/i?' 1 ' 

and the corresponding timescale for damping due to bulk 
viscosity 
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with r sa 5 to simulate an almost uniform density. 

The growth time due to emission of gravitational waves 
(r gw = —2EJ (dE I dt)\ gw ) is to lowest order in the angular 
velocity, fi, given by 
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If the timescale for gravitational wave emission is short 
compared to the damping timescales t sv and Tb v , the star 
will spin down. 

To find the critical angular velocity for a given stellar 
model as a function of temperature one solves the equa- 
tion 
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In the case of hot neutron stars, the / = 2, m = 2 r-mode 
instability was found to be decisive in the sense that it 
corresponds to the longest viscous timescales. 

Introducing R w = i?/10km, Mi. 5 = M/1.5M©, where 
Mq is the solar mass, and 5^3 = fi/10 3 s — 1 , one gets 
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For strange quark matter the viscous timescales are 



r sv (SQM) = 1.01 x I0 8 s (as /0.1) 5 / 3 MZl /9 Rll /3 T^ s , 

(9) 



T bv (SQM) = 5.75 x 10~ 2 s', 
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Here mioo = m s /100MeV, and for simplicity only the In- 
dependent term, not the Aw/w -term in Eq. (||) has been 
included. 

For comparison the shear viscosity in ordinary neutron 
star matter is dominated by electron-electron scattering 
for Tg < 1. The corresponding timescale is [|J 
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r sv (NS) = 3.4 x 10 7 sAfi:^ lf . 
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If the proton fraction in the neutron star is below g , the 
modified URCA-process is the dominant weak interac- 
tion, giving a bulk viscosity damping time H 



Tbv(NSmU) = 4.33 x 10 9 s njM££iq a Tg 
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If a larger proton fraction is available, as happens in some 
equations of state, the faster direct URCA-process is ac- 
tive fl7j], and the bulk viscosity timescale changes to 



T bv (NSdU) = 227s(<jx ir 1/3 ^f 1 2/ 5 3 Tr 4 



(13) 



where qx i » 1 (c.f. |Tq]). 

Inserting these timescales in Eq. (^) one can solve for a 
critical angular velocity, Q c , for the onset of the r-mode 
instability. Figure 1 shows the results as functions of 
temperature for stellar models with M1.5 = Rio = 1. fi c 
is expressed in units of the maximum (Kepler) angular 

,l/2 p -3/2 
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velocity, Q, K » 0.67v^FGp w Slabs' 1 Ml',? R~i 



For neutron stars with proton fraction below ^ the 
transition between shear- and bulk viscosity damping is 
seen to take place at T w 10 9 K, whereas direct URCA- 
processes reduces the transition temperature to 10 8 K. 
The time required to cool the interior of a neutron star 
to 10 9 K by the modified URCA mechanism is less than 
a year, and via the direct URCA mechanism, the cooling 
timescale is roughly t coo \ rs 20sT 9 -4 (see J01 for a re- 
view of neutron star cooling), so in any case, the neutron 
star will be subject to significant r-mode damping of the 
rotation rate within a matter of days to months. 

In contrast, the transition between shear- and bulk 
viscosity damping in a strange star only takes place af- 
ter the interior of the star has cooled to around 4.6 x 
10 6 K(a s /0.1) 10/19 m ro 1 8/19 , and with a reduction in fl c 
to between 0.26 and 0.32 of the Kepler frequency for a 
strange quark mass of 100 and 200 MeV respectively, 
corresponding to rotation periods of 3 milliseconds or 
faster, rather than the much stronger reductions expe- 
rienced by a neutron star. Furthermore, the time re- 
quired for a strange star to cool to the temperatures 
where the r-mode instability rotation damping sets in 
is much longer than the few months in the case of a neu- 
tron star. A typical cooling timescale for a strange star 
is icooi ~ 10 -4 yearsXg -4 , so that cooling to 4.6 x 10 6 K 
would take 2 x 10 5 years (and if the electron fraction in 
the star is very low, cooling may have to proceed on a 
much longer timescale, the quark matter analog of the 
modified URCA process; see po) for a recent study of 
strange star cooling). 

If the initial rotation rate of the strange star is very 
rapid, say half the Kepler frequency, the r-mode insta- 
bility will set in after a few years for low strange quark 
mass, but increasing the rotation period even above 2 ms 
will take more than 10 3 years (and to reach the ultimate 
limit of 3 ms takes at least a hundred times longer). 



These results were derived keeping only the temper- 
ature dependent term in the bulk viscosity of strange 
quark matter, Eq. (0) p3]. The actual viscosity is even 
higher due to the presence of the nonlinear term depend- 
ing on the amplitude of the density perturbation. Thus 
the bulk viscosity will dominate at even lower tempera- 
tures than found here, but a detailed numerical study of 
the eigenmodes in a stellar model is required to include 
the nonlinear term properly. 

It therefore seems safe to conclude, that strange stars 
differ significantly from ordinary neutron stars in terms 
of their (lack of) sensitivity to the spin-down due to 
the r-mode instability. A young, rapidly rotating pul- 
sar with period below 5-10 milliseconds could therefore 
very well be a strange star, if strange quark matter is 
absolutely stable, or a neutron star with a significant 
core of strange quark matter, if strange matter is only 
metastable (a so-called hybrid star) 22 1. At present no 
such young millisecond pulsars (i.e. young pulsars with 
periods below 10 milliseconds) are known, which may ei- 
ther indicate, that pulsars are neutron stars rather than 
quark stars, or perhaps that the stars are all formed with 
low angular momenta for some other reason. Old neu- 
tron stars being spun up by angular momentum transfer 
in binaries are not subject to the r-mode instability, so 
the presently known millisecond pulsars could be either 
neutron or strange stars, whereas an old submillisecond 
pulsar because of the extremely high bulk viscosity of 
strange quark matter might be a strange star but prob- 
ably not a neutron star jl4],[l6 10 1. 

More detailed calculations will be necessary to set the 
exact limits for the rotation frequency of neutron stars, 
hybrid stars and quark stars as a function of age and 
cooling history, but the discovery of even a single young 
pulsar in the millisecond regime would be most exciting 
for the strange quark matter hypothesis. 



ACKNOWLEDGMENTS 

This work was supported in part by the Theoretical 
Astrophysics Center under the Danish National Research 
Foundation. 



[11 N. Anderss on, Ap. J. (1998) in press (preprint gr-| 
qc/9706075| ). 



[2] 
[3] 
[4] 



J. L. Frie dman and S. M . Morsink, Ap. J. (1998) in press 
(preprint ^r-qc/970607^ ). 

L. Lindblom, B. J. Owen, and S. M. Morsink, Phys. Rev. 
Lett. (1998) in press (preprint |gr-qc/9803053l ). 
N. Andersson, K. D. Kokkotas, and B. F. Schutz, preprint 
astro-p h/9805225| . 



3 



[5] 

[6] 
[7] 



K. D. Kokkotas and N. St ergioulas, Astron. Astrophys. 
(1998) submitted (preprint |astro-ph/9805297|) . 
B. J. Owen et al, preprint gr-qc/9804044 . 
For a Newtonian, uniformly rotating, uniform density 
star the maximal angular velocity is Qk ~ 0.67^-kGp ~ 



[8] 
[9] 

[10] 
[11] 



[12] 



[13] 
[14] 
[15] 

[16] 
[17] 

[18] 
[19] 
[20] 

[21] 



[22] 



8185s" 1 (Af/1.5Af Q ) 1/2 ( fi/10km 



-3/2 



N. Stergioulas, preprint gr-qc/9805012. 

H. Heiselberg and C. J. Pethick, Phys. Rev. D 48, 2916 

(1993). 

J. Madsen, Phys. Rev. D 46, 3290 (1992). 

E. Witten, Phys. Rev. D 30, 272 (1984); P. Haensel, J. L. 

Zdunik, and R. Schaeffer, Astron. Astrophys. 160, 121 

(1986); C. Alcock, E. Farhi, and A. Olinto, Ap. J. 310, 

261 (1986); N. K. Glendenning, Compact Stars, Springer 

(1997). 

The rate for the reaction u + d <-> s + u was re-calculated 
in J. Madsen, Phys. Rev. D 47, 325 (1993), and in H. 
Heiselberg, Phys. Scr. 46, 485 (1992). Earlier calcula- 
tions, including H. Heiselberg, J. Madsen, and K. Riis- 
ager, Phys. Scr. 34, 556 (1986), were wrong. 
Q. D. Wang and T. Lu, Phys. Lett. 148B, 211 (1984). 
R. F. Sawyer, Phys. Lett. 233B, 412 (1989). 
C. Cutler, L. Lindblom, and R. J. Splinter, Ap. J. 363, 
603 (1990). 

M. Colpi and J. C. Miller, Ap. J. 388, 513 (1992). 

J. M. Lattimer, C. J. Pethick, M. Prakash, and P. 

Haensel, Phys. Rev. Lett. 66, 2701 (1991). 

J. L. Zdunik, Astron. Astrophys. 308, 828 (1996). 

C. J. Pethick, Rev. Mod. Phys. 64, 1133 (1992). 

C. Schaab, B. Hermann, F. Weber, and M. K. Weigel, 

Ap. J. 480, Llll (1997). 

Z. Dai and T. Lu, Z. Phys. A 355, 415 (1996) (see also A. 
Goyal et al, Z. Phys. A 349, 93 (1994)) have estimated 
the influence of as on the bulk viscosity. The corrections 
found are too small to influence the conclusions reached 
above for realistic values of as- 

Another interesting possibility would be that some other 
exotic phase in the neutron star led to a dramatic in- 
crease in the viscosity. The direct URCA process is an 
example of a process that does increase the bulk viscos- 
ity, but in a less dramatic way. Nuclear matter with high 
strangeness content or even more exotic phases might ex- 
perience something similar. However, since such phases 
are typically concentrated near the center of the star, 
nothing as dramatic as seen for a strange star with high 
viscosity quark phase almost to the surface should be 
expected. 
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FIG. 1. Critical angular frequencies in units of the Ke- 
pler frequency as functions of temperature derived from Eq. 
([?]). Dotted curve corresponds to the shear viscosity for neu- 
tron stars, dashed curve to shear viscosity for strange stars. 
Dashed-triple-dot curve includes bulk viscosity from modified 
URCA processes in neutron stars; dash-dot curve bulk vis- 
cosity from direct URCA processes in neutron stars, and full 
curves are based on bulk viscosity of strange quark matter 
for m s — lOOMeV (lower curve), and m s = 200MeV (upper 
curve). The r-mode instability is active when simultaneously 
above the relevant curves for shear and bulk viscous damping. 
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